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ABSTRACT 

Context. Most of the optically classified low ionisation narrow emission-line regions (LINERs) nuclei host an active galactic nuclei 
(AGN). However, how they fit into the unified model (UM) of AGN is still an open question. 

Aims. The aims of this work are to study at mid-infrared (mid-IR) (1) the Compton-thick nature of LINERs (i.e hydrogen column 
densities of Nr > L5 x lO^^^cm"^); and (2) the disappearance of the dusty torus in LINERs predicted from theoretical arguments. 
Methods. We have compiled all the available low spectral resolution mid-IR spectra of LINERs from the InfraRed Spectrograph (IRS) 
onboard Spitzer. The sample contains 40 LINERs. We have complemented the LINER sample with SpitzerjlKS spectra of PG QSOs, 
Type-1 Seyferts (Sis), Type-2 Seyferts (S2s), and Starburst (SBs) nuclei. We have studied the AGN versus the starburst content in our 
sample using different indicators: the equivalent width (EW) of the polycyclic aromatic hydrocarbon (PAH) at 6.2pm, the strength 
of the silicate feature at 9.1pm, and the steepness of the mid-IR spectra. We have classified the spectra as SB-dominated and AGN- 
dominated, according to these diagnostics. We have compared the average mid-IR spectra of the various classes. Moreover, we have 
studied the correlation between the 12pm luminosity, yLy(12yL/m), and the 2-10 keV energy band X-ray luminosity. Lx(2 - 10 keV). 
Results. In 25 out of the 40 LINERs (i.e., 62.5%) the mid-IR spectra are not SB-dominated, similar to the comparison S2 sam¬ 
ple (67.7%). The average spectra of both SB-dominated LINERs and S2s are very similar to the average spectrum of the SB class. 
The average spectrum of AGN-dominated LINERs is different from the average spectra of the other optical classes, showing a 
rather fiat spectrum at 6 - 2^pm. We have found that the average spectrum of AGN-dominated LINERs with X-ray luminosi¬ 
ties Lx(2 - lOkeV) > 10"^^erg/s is similar to the average mid-IR spectrum of AGN-dominated S2s. However, faint LINERs (i.e. 
Lx(2 - 10 keV) < lO'^^erg/s) show fiat spectra different from any of the other optical classes. The correlation between vLy(l2pm) and 
Lx(2 - 10 keV) for AGN nicely extends toward low luminosities only if SB-dominated LINERs are excluded and the 2-10 keV band 
X-ray luminosity is corrected in Compton-thick LINER candidates. 

Conclusions. We have found that LINERs proposed as Compton-thick candidates at X-ray wavelengths may be confirmed according 
to the X-ray to mid-IR luminosity relation. We show evidence in favour of the dusty-torus disappearance when their bolometric 
luminosity is below Lboi - lO'^^erg/s. We suggest that the dominant emission at mid-IR of faint LINERs might be a combination of an 
elliptical galaxy host (characterised by the lack of gas), a starburst, a jet, and/or ADAF emission. Alternatively, the mid-IR emission 
of some of these faint LINERs could be a combination of elliptical galaxy plus carbon-rich planetary nebulae. In order to reconcile the 
Compton-thick nature of a large fraction of LINERs with the lack of dusty-torus signatures, we suggest that the material producing 
the Compton-thick X-ray obscuration is free of dust. 
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1. Introduction 

The emission in active galactic nuclei (AGN) is powered by ac¬ 
cretion onto a supermassive black hole (SMBH). AGN are tra¬ 
ditionally divided into two main classes based on the presence 
(Type-1) or not (Type-2) of broad permitted lines (EWHM>2000 
km s“^) in the optical spectrum. The so-called unification model 


^ Juan de la Cierva Fellow and Tenure track at CRyA (e-mail: 
o.gonzalez@crya.unam.es) 


(UM) proposes that both types of AGN are essentially the same 
objects viewed at different angles ( |Antonuca]|I993t |Urry & 


|Padovani||I995| ). An optically thick dusty torus surrounding the 
central source would be responsible for blocking the region 
where these broad emission lines are produced (the broad line 
region, BLR) in Type-2 Seyferts. The torus must not be spheri¬ 
cally symmetric, in order to obscure the BLR, while allowing at 
the same time the emission coming from region producing the 
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permitted narrow lines (known as narrow-line region, NLR) to 
reach us from the same line of sight (LOS). 

Low-ionisation nuclear emission-line regions (LINERs), first 
classified by |Heckman| ( |1980' I, are the dominant population of 
AGN in the local Universe ( |Ho et ar]|1997| ). However, they re¬ 
main as one of the most captivating subsets of nuclear classes 
because their main physical mechanism is still unknown. The 
nature of LINERs was initially sustained in their optical spec¬ 
trum, which can be reproduced with a variety of different phys¬ 
ical processes (e.g. photoionisation from hot stars, non-thermal 
photoionisation, shocks, post-main sequence stars, or AGN, 


pita & Sutherland] 1995 |Heckman|1980t|Eerland & Netzer|1983 
Veilleux & Qsterbrock 1987t|Stasinska et al.|2008t Singh et al 


([Gonzalez-Martm et al.|2006 

2009 Al Dudik et al. 20091 Younes 

et al.l|2011l 

Asmus et al.[[20] 

[11 [Mason et al.[[2012l [Gonzalez- 

Martin et al. 

2014[) ^ This is confirmed also from X-ray and UV 

variability studies ([Maoz et al. 20051 Hernandez-Garcia et al.[ 


obscuration is responsible for their differ ences (e.g. [Dudik et 


al. 


al. 


|2009t [Gonzalez-Martm et aL]|2009B| ). |Gonzalez-Martm et 


as a tracer of Compton-thickness (i.e. Nh > 1.5 x 10"^"^ cm ^), 


Gonzalez-Martm et^ ( |2009B] ) found that up to 53% of the LIN- 


does not suffer such a large extinction. Eurthermore, the dust 
that absorbs the shorter wavelength emission reradiates in the 
mid-IR and correspondingly produces a substantial fraction of 
the bolometric flux of the object. Dissecting the detailed mid-IR 
spectra of AGN can reveal the properties of the dust in the nu¬ 
clear region (e.g. |Mendoza-Castrej6n et al.|2015| ). Subarcsecond 
resolution studies claim a tight correlation between rest frame lu¬ 
minosities yLy(12yum) and Lx(2 - 10 keV) for Type-1 and Type- 


2 Seyferts ([Horst et al. [[20091 [Gandhi al.[[2009l 

Asmus et al. 

[201 ![ [Masegosa et al.[20 13 [[Gonzalez-Martm et al. 

201 3[). Imag- 


2013| ). In fact, many authors have suggested that LINERs are 
an he terogeneous class ( [Satyapal et al.|2005} [Dudik et al.|2005| 
2009| ). Indeed the emission mechanism dominating their optical 
spectrum is still under debate. [Singh et ah] p013| ) show that an 
AGN alone cannot explain the radial profile of the surface bright¬ 
ness Ha emission line; a contribution of extended emission post- 
main sequence stars is needed at optical frequencies. The reason 
why this contribution can be seen in LINERs might be the intrin¬ 
sic weakness of the AGN, which would outshine these signatures 
in more powerful AGN. In line with that, Gonzalez-Martm et al.[ 
( 2014[ ) showed that the host-galaxy contributes in a large frac¬ 
tion in most of the LINERs even at X-rays. Thus, the analysis of 
LINERs could be key to study the interplay between the AGN 
and the host galaxy. 

Nowadays we know that around 75-90% of LINERs 
show evidence of AGN using multi-wavelength information 


ing data at mid-IR wavelengths have shown that some Compton- 
thick LINER candidates might fa ll off this relation, with Li 2 ^m 
larger than expected for their Lx ( [Mason et al.|2012[ [Masegosa 


other AGN? How do they fit into the UM of AGN? Some re¬ 
sults suggest that they might constitute a class of AGN with 
a different accretion mode (e.g. [Younes et al.[[2Qll| [Nemmen| 


et al. 2014) while some other authors have argued that large 


_ ( [2009A[ ) found that the hydrogen column density, Nh, in 

LINERs range from the galactic value up to Nh - 10^"^ cm“^. 
This is fully consistent with the Nh values reported for Seyfert 
galaxies (e.g. [Panessa et al.|[2006[ [Bianchi et al.[[2012l [Mar-[ 
|inucci et al.[[2012 ). However, using the ratio between the lu¬ 
minosity of the [OIII]/15007A emission line and the intrinsic 
hard (2-10 keV) X-ray luminosity, L([OIII])/Lx(2 - 10 keV), 


ERs in their sample are Compton-thick candidates. This percent¬ 
age is two times higher than that reported for Type-2 Seyferts 
([Maiolino et al.[[T99^ Bassani et al. [1 1999| [Panessa et al.|2006t 


^ Note that most of these studies select their sources using archival X- 
rays observations so they might be biased due to the complex selective 
effect this introduces in the sample. 


[et al.[[20l3] ). The confirmation of such a trend for a large sam¬ 
ple might confirm the Compton-thick nature of a large portion 
of LINERs. 

Erom the theoretical point of view, [Elitzur & Shlosman 


( [2006| ) showed that the torus might disappear when the bolomet¬ 
ric luminosity decreases below Lboi ~ lO'^^erg/s because the ac¬ 
cretion onto the SMBH cannot longer sustain the required cloud 
outflow rate. Thus, the low bolometric luminosity of LINERs 
makes them key to probe this theoretical prediction. [Maoz et al. 


( [2005 1 ) showed that the fraction of variable Type-1 and Type-2 
LINERs at UV is the same. This favours the lack of a dusty torus 
obscuring the central AGN in Type-2 LINERs since the central 
source is the responsible for such variability pattern. The mid- 
IR emission shows clear signatures of the dusty torus. In fact. 


explaining the mid-IR emission of Seyfert galaxies (e.g. 

, Ramos 

Almeida et al.|2009 201 

Alonso-Herrero et al.|2011l I 

ionig et 

al.|2010[). Thus, the mid-IR spectra of LINERs can give impor- 


tant clues on the existence of the dusty torus for low bolometric 
luminosities. However, other contributors like jet emission (e.g. 
NGC 1052, [Eemandez-On tiveros e t al.|[2Q12[ ) or optically thin 


dust (e.g. NGC 3998, [Mason et al.|[2Q13| ) can also contribute to 


the mid-IR emission. 

The purpose of this paper is to study (1) the Compton-thick 
nature of LINERs and (2) the plausible disappearance of the 
torus. We present the mid-IR Spitzer/IRS spectra of 40 LINERs. 
We compare them with mid-IR Spitzer/IRS spectra of Seyferts, 
PG QSOs, and Starbursts (see Sectionj^ The data reduction and 
measurements are described in Section]^ The relatively low spa¬ 
tial resolution of Spitzer/IRS spectra makes these spectra being 
contaminated from the host galaxy emission, which is partic¬ 
ularly relevant for low-luminosity AGN (LLAGN) as LINERs. 
Section[4]describes a method which is able to select those mid-IR 
spectra with a negligible starburst contribution. Section stud¬ 
ies the average spectrum of LINERs and compare them with that 
of Seyferts, PGQSOs, and Starbursts. Sectionj^ shows the anal¬ 
ysis of the correlation between Lx(2 - 10 keV) and Lujjm foi* 
LINERs. In Section [7] we discuss the implications of the main 
results. The conclusions of this paper are summarised in Section 


[Cappi et al.|20Q6| ). [Dudik et al.[ ( |2009| ) studied the emission lines 

in 67 high-resolution Spitzer/IRS spectra of LINERs and found 
that the central power source in a large percentage of LINERs is 
highly obscured at optical frequencies, consistent with the X-ray 
results. 

Obscuring dust hampers the studies of the optical to soft X- 
ray emission coming from both young hot-stars and the accretion [^1 ^1- 
disk. However, emission at mid-infrared (mid-IR) wavelengths 


2. Sample 

Our initial sample of LINERs comes from the catalog of 
LINER s obser ved at X-rays published by [Gonzalez-Martm 
( 2009A[ ). This guarantees that all the LINERs have 


Lx(2 - 10 keV) measurements, what is crucial for our purposes. 
However, we must be aware that this sample does not constitute 
a complete sample. Eor the present analysis we have used two 
databases to obtain the Spitzer data for an additional sample of 
LINERs. The first one is the Cornell atlas of Spitzer/IRS spectra 
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(CASSIsj^. CASSIS provides low-resolution spectra (R~60-I27 
over 5.2yum to 38yum) with the IRS instruments in the stare mode 
( [Lebouteiller et al.||20II| ). The second database is the Spitzer 
infrared nearby galaxy survey (SINGS [Kennicutt et al.|[200^ . 
SINGS is a Legacy programme of imaging and spectroscopic 
data for 75 nearby galaxies. As part of the Legacy programme 
all the one-dimensional nuclear spectra have been archived in 
the infrared science archive (IRSv^. This provides uniform 5- 
30/im spectra in the spectral mapping mode. Most galaxies in the 
SINGS sample have also been observe d with Chandra an d their 
main X-ray properties are described by [Grier et al.| ( |20lT] ). 

The CASSIS atlas conta ins 27 LINERs published by 
[Gonzalez-Martm et~aL] ( |2009A|). We have also adde d NGC 3079, 
that was not analysed by Gonzalez-Martm et al.[ ([2009A[) but 
it was included as a Compton-thick LINER by [Goulding et al 


([20I2[). Grier et aL]([20II[) inc luded 20 LINERs (8 included in 
[Gonzdlez-Martm et al.|2009A| ). 

The final sample of LINER in this paper contains 40 
Spitzer spectra. Among them eight have been optically clas¬ 
sified as Type-1.9 LINERs (LINER l) and 32 as Type-2 LIN¬ 
ERs (LINER2) by [Ho et ( [1997[ ). Three of them are known 
Compton-thick, 21 are Compton-thin, and 16 were classified as 
Compton-thick candidates by [Gonzalez-Martm et al.[ ([2009B[). 

[Dudik et ( [2009[ ) reported mid-IR spectra of 67 LINERs 
(13 objects in common with our sample). However, they used 
the high-resolution modes of Spitzer because their work was 
focused in the fine structure mid-IR emission lines. Thus, their 
results are not directly comparable with ours. [Sturm et al.[p006[ ) 
reported a mid-IR study of 33 LINERs. They selected their sam¬ 
ple on the basis of IR luminosity while our sample is conformed 
by objects with measured X-ray luminosities. As a consequence, 
only NGC 4486 is common with our sample. 


2.1. Comparison samples 

To study the possible contribution of star-formation or AGN 
emission, we have selected starbursts, Seyferts and PG QSOs to 
be able to compare their mid-IR spectra to those of LINERs. 
Note that the sample is not complete in any sense but it allows 
us to have a representative set of objects for each category to 
compare LINERs with them: 


- Seyferts. All the Type-1 and Type-2 sources included in[Shi| 


et al. 

(2006), in 

Goulc 

ing et al. 

samp 

e. In tota 


( [2012[ ), and those included in the SINGS 
it contains 42 Seyferts. Among them 31 
are Type-2 Seyferts (S2, including 19 Common-thick and 12 
Compton-thin) and 11 Type-1 Seyferts (Sl^ 

Palomar Green QSOs (PGQSOs). Thi s sample includes all 
the PG QSOs in the sample defined by [Veilleux et al. 


with Spitzer/IRS spectra in CASSIS and redshifts z <0.2^1 
This PGQSO sample includes 26 sources. 

Starbursts. Th i s samp le is t aken from[Ranalli et al.[ ( [2003[ ), 
[Brandi et al.[ ( [2006[ ), and [Grier et al.[ ( [2011[ ). The star- 
burst sample contains 21 sources. Note that among them 
NGC 3367 was classified as a Seyfert by [Veron-Cetty &[ 


^ http://cassis.astro.cornell.edu/atlas/ _ 

^ http://irsa.ipac.caltech.edu 

Note that the SI sample contains objects classihed as I’ype 1, I’ype 

1.2, Type 1.5, Type 1.8, a nd Type 1.9 Seyferts. The S2 sample in clud es 
only purely Type-2 Seyferts. _ 

^ Ihe redshift limit has been chosen to be able to obtain rest frame 

30yum luminosities, required for our analysis. 


Veron[([2006[) although it was classified as a starburst by[Ho[ 
et al.[ ( [1997[ ). 


Note that for all the classes we have only included spectra ob¬ 
served with both the short-low (SL) and long-low (LL) modules 
to guarantee the full Spitzer/IRS coverage (at least ~ 5 - 30yum). 
Moreover, ultra-luminous infrared galaxies (ULIRGs) have been 
excluded from the analysis because they might have a controver- 
sial source of emissio n at mid-IR pmanishi et al.|2007] [Alonso- [ 
[Herrero et al.[[201^ . All together these samples comprise 89 
sources (129 nuclei including LINERs). 


3. Data processing and analysis 

CASSIS and SINGS provide fiux and wavelength calibrated 
spectra. However, the observations using data from both the SL 
and LL spectral modules suffer from mismatches due to tele¬ 
scope pointing inaccuracies or due to different spatial resolution 
of the IRS orders. This is not corrected in the final products given 
by CASSIS and SINGS. We therefore scaled each spectra to the 
immediate prior (in wavelength range) to overcome such effects. 
Thus, our fiux level is scaled to the level of the shortest wave¬ 
lengths, which is the order with the highest spatial resolution. 
This guarantees that the fiux level is scaled to the best spatial 
resolution that Spitzer can provide. Moreover, the spectra are 
shifted to rest-frame according to the redshift of the objects. 

Eor each object we have measured the \2pm and 30yum lu¬ 
minosities using the Spitzer/^RS spectra. Errors have been esti¬ 
mated assuming 15% flux-calibration uncertainties, which fully 
dominate other source of errors (e.g. [Gonzalez-Martm et al.[ 
[2013t [Ramos Almeida et al.|2011[ ). We have also measured the 
fluxes and equivalent width (EW) of the polycyclic aromatic hy¬ 
drocarbon (PAH) features at 6.2 and 11.3 pm. The EW of the 
PAH features were measured by integrating the emission over 
the continuum in a wavelength range of 5.9-6.5pm and 11.0- 
ll.6pm for the 6.2 and 11.3yum PAH emission features, respec¬ 
tively. The continuum was estimated through a linear fit to the 
5.5-5.9yum (10.7-1 l.Oyum) and 6.5-6.1pm (11.6-11.9yum) for the 


6.2pm (11.3yum) PAH feature (see e.g. [Gonzalez-Martm et al.[ 
[2013[ ). We have also computed the strength of the silicate emis¬ 
sion/absorption feature at 9.7pm through the apparent depth at 
9.1pm, T 9 . 7 ^ni (e.g. Shi et al.|2006[ Levenson et al.|2007 ): 


'r9.1iim ~ cont,9.1 jim/7^9Jinn) ( 1 ) 

where p 9 . 7 ;^m and Pcont, 9 . 7 /im are the fluxes of the spectra around 
9.1pm and its expected continuum, respectively. Note that the 
apparent depth at 9.1pm T 9 . 7 ^^ is positive for absorption silicate 
features and negative for emission features. 

Due to the complexity of Spitzer spectra, we have used PAH- 
Pnj^to obtain T 9 . 7 ^^ except when emission silicate features are 
detected (see below). PAHEIT is an IDL tool for decompos¬ 
ing Spitzer/IRS spectra of PAH emission sources, with a spe¬ 
cial emphasis on the careful recovery of ambiguous silicate ab¬ 
sorption, and weak, blended dust emission features ( [Smith et 
[2007[ ). PAHEIT is primarily designed for its use with the full 5- 
35pm Spitzer/IRS low-resolution spectra. However, PAHEIT is 
not able to treat or recover silicate emission features expected 
to occur in Type-1 AGN, giving T 9 . 7 ^^ = 0. In these cases we 
have computed T 9 . 7 ^^ by fitting the 9-14yum Spitzer/IRS spec¬ 
tra to a Gaussian profile. This is a general treatment to estimate 

^ http ://tir. astro.utoledo.edu/j dsmith/research/pahfit.php 
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Fig. 1. Histograms of the EW of the PAH at 6.2 pm (left) and at 11.3 pm (right) for PGQSO (top-panel, orange-broad bars), SI (top-panel, 
red-narrow bars), SB (middle-panel, green-broad bars), S2 (middle-panel, yellow-narrow bars), LINERl (bottom-panel, purple-narrow bars), and 
LINER2 (bottom-panel, pink-broad bars). The median values and 25th-75th percentile range for each class of object are shown with large crosses 
(with the same color-code than the histogram) with continuous lines for PG QSO, SB, and LINER2 and with dashed lines for SI, S2, and LINERl. 
The vertical dot-dashed line shows the limit chosen to divide the sample into strong- and weak-PAHs (see text). 


that has been proven to be a good approximation when 
compared with PAHFIT (see |Gonzalez-Martm et al.|2013| ). 

Together with the Spitzer/IRS spectra of the sample we also 
compiled the 12pm luminosities obtained with ground based 
telescopes. These measurements have the advantage to better 
isolate the nuclear source because they come from images with 


tribution might be very strong at mid-IR wavelengths, dominat¬ 
ing the entire emission ( Mason et al.|2012} . 


Over the last decade several diagnostics have been proposed 
to quantify the contribution of star formation and AGN activity 
to the infrared luminosity. These diagnostics are based on the 
mid-IR continuum slope, the EW of the PAH features, the ra¬ 
tio of [NeV] (or [OIV]) over [Nell], and the EW of the PAH at 


objects. 

We have also included four LINERs observed with 
GTC/CanariCam as proprietary data (program IDs GTC42-12B 
and GTC35-13A). We have reduced them uniformly with the 
RedCan package ( |Gonzalez-Martm et al.|2013| ). We present here 
their luminosities (marked with asterisks in Col. 5 in Table [TJ 
whereas the full imaging analysis will be presented in a forth¬ 
coming paper. All together we have ground-based measurements 
for 61 out of the 129 sources included in this paper. 

Table contains the X-ray luminosities, mid-IR measure¬ 
ments for the Spitzer/IRS spectra, and the 12pm luminosities 
from Spitzer/IRS and ground-based telescopes when available. 


nearby galaxies). Most of these measurements come from the 

6.2 or 11.2pm versus the 9.1pm optical depth ([Genzel et 

catalog of sub-arcsecond mid-IR images of AGN reported by 

al.jj 19981 [Lutz et al.jj 19981 [Dale et al.||2006t [Sturm et al.||2006t 

lAsmus et al.|(|20T3li. It includes Subaru/COMICS, VLTA^ISIR, 

Spoon et al. 

|2007| Baum et al.|2010| Heman-Caballero & Hatz-] 

Gemini/MICHELLE and Gemini/T-ReCS mid-IR data of 253 

iminaogloul 

201 1[). In this section we use several diagnostics to 


separate which Spitzer/IRS spectra are strongly contaminated by 
non-AGN emission. Note, however, that the starburst contribu¬ 
tion to the mid-IR spectra does not exclude the presence of an 
AGN. This is an attempt to determine whether the AGN is dom¬ 
inating or not the mid-IR spectrum. We have excluded from the 
analysis the diagnostics based on fine structure emission lines 
because they are blended with other emission lines at the spectral 
resolution of these dat aset. For more detail s in these diagnostics 
we refer the reader to [Dudik et al.| (2009), where they studied 
these emission lines for a large sampL 

4.1. PAH features 


e of LINERs. 


4. AGN versus starburst contents 

Spitzer has been used to study the largest samples of AGN ever 
analysed at mid-IR (e.g. |Shi et al.|2006[|Deo et al.|2007| ). How¬ 
ever, a disadvantage of these data is their relatively low spatial 
resolution. This makes Spitzer spectra to be often contaminated 
by the host galaxy. It is expected to be particularly relevant for 
LINERs where the AGN is faint. In this case the non-AGN con¬ 


The star formation activity correlates with the PAH strength, 
where starburst-dominated galaxies are then expected to show 
strong PAH features (e.g. [Peeters et al.||2004] ). This well stab- 
lished correlation has led to the use of PAH strength as a tracer 
of star formation (e.g. [Esquej et al.||2014|). PAHs might be de- 
stroyed due to the presence of an AGN ( Genzel et al.|1998[|Wu 


|et al.|2009] ). This is particularly relevant for the PAH feature at 
6.2pm that is produced by grains with smaller sizes and, there¬ 
fore, their destruction near the AGN is more efficient ( [Diamond^ 
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Fig. 2. The silicate apparent depth at 9.7 //m, versus the EW of PAH at 6.2 //m, EW(PAH 6.2 //m), for PGQSOs, Sis, and SBs (left) and for 
LINERs and S2s, (right). LINERls, LINER2s, Sis, S2s, PG-QSOs, and starbursts are shown with purple and pink circles, red and yellow squares, 
orange stars, and green up-side down triangles, respectively. White circles and white stars mark known Compton-thick sources and Compton-thick 
candidates, respectively. Note that the error bars in T|^are always within the size of the symbol. In the right pane l we also include PG QSOs, Sis, 
and SBs with white symbols for comparison purposes. Green-dotted lines indicate the diagonal branch found by |Spoon et al.|p007|) for ULIRG s 
and Starbursts. The shadowed (grey) area shows the range of values for rjgjjthat could be explained with Clumpy models { [Nenkova et al.|2008| ). 
We also show with orange and yellow li nes the expected range of values for rj^in face-on AGN (assuming i = (J^ and in edge-on AGN (assuming 
i = 9C(^ using the models described by |Nenkova et al.| ( |2008| ) (see text). The red-short vertical line shows EW(PAH = 0.228yum which 

divides into weak and strong PAHs. 


|Stanic & Rieke||2012| ). However, the ll.Syum might not be sup¬ 
pressed by the AGN (at distances as close as 10 pc) but diluted 
when the AGN continuum becomes dominant ( |Alonso-Herrero| 
|et al.||2014l [Ramos Almeida et~aL||2Q14| ). Note that our aim is 
to select those spectra where the host galaxy contribution due 
to star-formation is not dominating the mid-IR spectrum. The 
PAH strength is a good tracer of star-formation occurring at far 
distances from the AGN, where this destruction/dilution of the 
PAH features is negligible. Supporting this, the PAH at 11.3//m 
was negligible in 18 out of the 20 AGN reported by [Gonzalez^ 
[Martin et al.j p013[ ) with high-spatial resolution spectra while 
their Spitzer spectra showed strong PAH features. Thus, these 
diagnostics are still useful in our analysis, irrespective of the di¬ 
lution or suppression of the PAH features near the AGN. 

Fig. shows the distributions of the EW(PAH) per class for 
the PAH features detected at 6.2/im (left) and \\3pm (right). 
Note that upper limits on the non detected lines are not included 
in these histograms. The median EW(PAH 62pm) for the SB 
class is significantly higher than that for the SI and PGQSO 
classes. Only one SI shows an EW(PAH 6.2yum) consistent 
with the SB class (NGC 5033). Moreover, only one SB shows 
a limit on the EW(PAH 62pm) consistent with Sis or PG QSOs 
(NGC 3184). The S2 class have objects with EW(PAH 62pm) 
overlapping with values found for Sis, PG QSOs, and SBs. LIN¬ 
ERs (both LINERl and LINER2), like S2s, also show EW(PAH 
6.2yum) spreading a large range of values. The histogram of 
EW(PAH \ \3pm) is similar to that of EW(PAH 62pm). How¬ 
ever, there is a larger overlapping between the distributions in 
EW(PAH 11.3yum) for SBs with Sis and PGQSOs than in the 


histogram of EW(PAH 6.2pm). Moreover, the S2 class overlaps 
with the SB class for a larger number of objects in EW(PAH 
11.3yum) than for the histogram of EW(PAH 6.2/im) (15 and 6 
objects, respectively). 

The PAH feature at 11.3/im might be strongly attenuated by 
the silicate absorption feature ( [Brandi et al.[[2006'| ). [Gonzalez^ 
[Martin et~ar] ( [20 13 [ ) estimated this attenuation to be up to -40% 
of the intrinsic PAH feature at 11.3/im for T 9 . 7 ^rn = 1- This per¬ 
centage can be higher for larger The EW(PAH 6.2/im) 

might be a better tracer of star-formation when the silicate at¬ 
tenuation is large because it is not embedded in the silicate ab¬ 
sorption feature. This is clearly seen in Fig.[^ where the overlap 
of the S2 class (expected to be more attenuated than Sis under 
the unified model) with the SBs is much higher for the EW(PAH 
11.3/im) than for the EW(PAH 6.2pm). We therefore have cho¬ 
sen EW(PAH 6.2/im) as a better tracer of the SB-dominance in 
our sample. 

We have defined a limit on EW(PAH 6.2/im) us¬ 
ing the mean value and the standard deviation over this 
mean value for objects classified as SBs as follows: 

< EW(PAH6.2/im) > -3 x o-(EW(PAH6.2/im)) = 0.233/im (i.e. 

< log(EW(PAH6.2/im)) >= -0.633). This ensures that 99.7% of 
SBs show EW(PAH 6.2/im) above this limit. Note here that the 
limit would be EW(PAH6.2/im) = 0.247 if the 99th percentile is 
used. This will give a slightly less restrictive limij^ Above (be¬ 
low) this value we classified the objects as strong-PAH (weak- 


^ Only two objects will be included in the weak-PAH class if the limit 
is set to EW(PAH6.2//m) = 0.247 compared to those obtained using 
EW(PAH6.2//m) > 0.233, namely NGC 7130 and NGC 3367. 
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PAH) objects. Six out of the 31 S2s are consistent with the 
strong-PAH category; ten out of the 40 LINERs are classified 
within the strong-PAH class. All of them are LINER2s except 
N GC1097. _ 

[Spoon et al.| ( |2007| ) presented a mid-IR diagnostic of the 
AG N/ULIRG content based on Tgj versus EW(P AH 6.2pm) (see 
also |Hernan-Caballero & Hatziminaoglou|2011| ). The advantage 
of this diagnostic is that it takes into account the effects of strong 
obscuration of the nuclear source. They showed that galaxies 
are systematically distributed along two different branches: ( 1 ) 
a horizontal line with T 9.7 < 1 of continuum AGN-dominated 
to PAH-dominated spectra and (2) a diagonal line going from 
deeply obscured (high T9.7 and low EW(PAH 6.2pm) to PAH- 
dominated spectra (low T9.7 and high EW(PAH 6.2pm). Seyferts 
and QSOs are found exclusively on the horizontal branch with 
T9.7 < 1. The large majority of LIRGs and ULIRGs in [Spoon ^ 
^ ( [2007 [ ) are located in the diagonal line. Starburst are placed 
at the end of the two branches, with large EW(PAH 6.2pm) and 
T9.7 < 1. They argued that these two branches reflect a funda¬ 
mental difference in the dust geometry in the two set of sources. 
The horizontal branch could have a clumpy structure while the 
diagonal might be smooth. 

Pig. shows T9.7 versus EW(PAH 6.2pm) for SBs, 
Sis, and PGQSOs in the left panel and LINERs and 
S2s in the right panel. PGQSOs and Sis (T9.7 < 1 and 
EW(PAH 6.2pm) < 0.228yum) are clearly distinguished from 
SBs (EW(PAH 6.2pm) > 0.233/tm). This result is fully consis¬ 
tent with that reported by [Spoon et al.[ ( [2007[ ). Our diagram 
shows very few nuclei with deep silicate features and weak PAH 
features. This was also found by [Spoon et al.[ ( [2007] ) with only 
eight over the 160 objects in their sample belonging to this cate¬ 
gory. We have tested the use of the EW(PAH 11.3//m) instead of 
EW(PAH 6.2pm) in this diagram, finding a similar result. This 
was already reported by |Heman-Caballero & Hatziminaoglou| 
( 2Q11[ ) in a large sample of Spitzer/IRS spectra. 

There is a maximum T9.7 expected under the predictions of 
the clumpy torus models for AGN. Larger values of T9.7 can 
be interpreted as significant contamination from the host galaxy 
( [Alonso-Herrero et al.[201 H[Gonzalez-Martm et al.[2013[ ). In or¬ 
der to investigate this issue, we have cormuted T9.7 using a set of 
models within the libraries of CLUMP^ These consist on a set 
of spectral energy distributi ons (SEDs) using the A GN clumpy 
torus emission described by [Nenkova et al. ([2008 ). The param¬ 
eter ranges chosen are those reported by [Gonzalez-Martm et al.[ 
( |2013[ ). We have downloaded the SEDs for a width of the toroidal 
distribution cr = 45°, a ratio between the outer and the inner ra¬ 
dius of the torus rout/^int = 200 , an exponential slope of the radial 
distribution of clouds q = 2 , an optical extinction of the clouds 
ranging ry = 5 - 150, and a number of clouds along the equator 
of the torus Nq = 2 - 20 clouds. The number of clouds along the 
LOS, N, depends on the inclination angle, /, as N = Noe“^ . 
We refer the reader to [Gonzalez-Martm et al. (|2013[) for details 
on the selection of these parameters and to [Nenkova et al.[ ( [2008[ ) 
for the details on the modelling. 

We have computed T9.7 for these SEDs using the same 
methodology as for the Spitzer/IRS spectra reported here. Eig. 

shows as a grey area the minimum and maximum T9.7 found 
using these models (-1.13 < T9.7 < 1.25). Thus, objects with 
T9.7 > 1.25 are not expected under any clumpy torus model. 
We use T9.7 > 1.25 to classify an object as deep-silicate. We 
also show the range of T9.7 expected for face-on AGN assum¬ 
ing i = 0° (-0.96 < T9.7 < 0.62) and for edge-on AGN assuming 


^ http://www.pa.uky.edu/clumpy/ 
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Fig. 3. Histograms of \og{yLy{20pm)/yLy(30pm)) for PGQSO (top- 
panel, orange-broad bars), SI (top-panel, red-narrow bars), SB (middle- 
panel, green-broad bars), S2 (middle-panel, yellow-narrow bars), 
LINER 1 (bottom-panel, purple-narrow bars), and LINER2 (bottom- 
panel, pink-broad bars). The mean values and one standard deviation 
over the mean for each class of objects are shown with large crosses 
(with the same color-code than the histogram), continuous lines for 
PG QSO, SB, and LINER2 while dashed lines for SI, S2, and LINERl. 
The grey area of the plot shows the range of values expected for AGN 
according to the models given by [Nenkova et al.[ ( |2008j ). The orange 
and yellow vertical lines show the same ranges but for inclination an¬ 
gles of i = 0° and i = 90°, assuming that these values are representative 
of face-on and a edge-on AGN (see text). 


i = 90° (-0.86 < T9.7 < 1.25). Sis are naturally explained within 
the range of values of T9.7 expected under the clumpy torus mod¬ 
els. S2s tend to show larger T9.7 than Sis. Only two S2s are out 
of the expected range with clumpy torus models. Only two SBs 
and two LINERs show T9.7 > 1.25. Eollowing this diagram, we 
have divided our sample into three categories: 

- Deep-Silicate: Strength of the silicate feature above the max¬ 
imum explained by clumpy models (i.e. T9.7 > 1 . 25 ) regard¬ 
less ofEW(PAH 6.2pm). 

- Strong-PAH: T9 7 < 1.25 and EW(PAH 6.2pm) > 0.233yum. 

- Weak-PAH: T9.7 < 1.25 and EW(PAH 6.2pm) < 0.233yum. 

This classification is included in Table [T] Starburst are lo¬ 
cated mostly in the region of strong-PAH. Two of them, though, 
populate the area of deep-silicates. PG QSOs and Sis are placed 
in the region of weak-PAH^ 

S2s tend to show larger T9.7 than Sis, as expected under the 
unified model of AGN. Similarly to Sis and PG QSOs, most S2s 
are in the region of weak-PAH. Only two of them are within 
the area of deep-silicates and five of them are within the area of 
strong-PAHs. Thus, according to this diagram, only seven out of 
the 31 S2s show signs of host-galaxy contamination at mid-IR. 

^ The only exception is NGC 5033 which is located in the region of the 
diagram of strong-PAHs. 
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Fig. 4. The 12yum flux obtained with ground-based telescopes (nuclear) versus 12yum flux obtained with the Spitzerj^R^ spectra. Note that both 
quantities are shown in log-scales. A cross is shown to illustrate the error bars in these measurements (see Section |^. The left panel shows the 
results for the entire sample of objects with ground-based telescopes, the middle panel shows the objects classifled as SB-dominated and the right 
panel shows only objects with Spitzerj^R.^ spectra classifled as AGN-dominated in this work (see text). The blue continuous line represents the 
one-to-one relation and the green dashed line the best linear relation. 


LINERs mostly populate the area of weak-PAHs. All the 
LINER Is but NGC 1097 are in this area of the diagram. Among 
LINER2s, two are in the region of deep-silicates and nine are in 
the strong-PAHs area. If LINERs in the weak-PAH area of the 
diagram are considered as AGN-dominated at mid-IR, then 30 
out of the 41 are AGN-dominated at mid-IR. 


4.2. Steepness of the mid-IR spectra 


The steepness of the mid-IR spectra characterises th e relative 
contribution of warm and cool dust to the mid-IR ( |Baum et] 
ar]|2010| ). Note that we refer here to the steepness at the spa¬ 
tial resolutions of the Spitzer/IRS spectra (i.e. kpc scales) while 
at smaller scales (below 100 pc scales obtained with ground- 
based instruments), the mid-IR emission is expected to be dom- 
inated from dust hea ted by the AGN ( |Honig et al.|2011t [Ramos | 
Almeida et al.|20lf] ). It has been proven as a good indicator of 
the starburst content given its correlation with EW(PAH 11 .3pm) 


(e.g.jWu et al.||2009l [We^man et al.|[2005| [Brandi et ar||2006 


LaMassa et al.]MT^ . This steepness has been defined in sev 

eral bands by different authors; e.g. 20-30yum ([Baum et al.poTO 


Weedman et al.|200^ or 1 5-30/im ( [Brandi et al.|2006[[Wu et al.[ 


2009 [Nardini et al.[[2008[ ). We present the 20 and 30pm lumi¬ 
nosities to produce an estimate of the steepness of the mid-IR 
spectra in our sample. 


Pig. shows the histogram of such steepness, expressed as 
log(yLy(20yum)/yLy(30yum)) (see also Table [^. PGQSOs and 
Sis tend to show larger values of log(yLy(20yum)/yLy(30yum)) 
than SBs. However, these two distributions overlap in the 
range -0.25 < log(yLy(20yum)/yLy(30yum)) < 0. S2s and LIN¬ 
ERs show a wide range of log(yLy(20yum)/yLy(30yum)), over¬ 
lapping with Sis, PGQSOs, and SBs, although Sis and 
PGQSOs distributions are skewed toward the larger values 
of log(yLy(20)Um)/yLy(30)Um)). Therefore, this ratio itself is 
not as good tracer of AGN dominance as it is EW(PAH 
6.2pm ) (see previous subsection). Although in theory it is a 
good tracer of the contribution of warm and cool dust to the 
mid-IR, in practice, some SB-dominated spectra can have a 
log(yLy(20yum)/yLy(30yum)) ratio consistent with those of AGN. 


Pollowing the same idea than in the previous subsection, we 
have computed log(yLy(20//m)/yLy(30yum)) of theoretical mod¬ 
els obtained with the Clumpy libraries. The vertical-thick orange 
and yellow areas (Pig. show the range of steepness expected 
from these models for face-on (assuming i = 0) and edge-on 
(assuming i = 90) torii, respectively (delimited also with verti¬ 
cal orange and yellow, respectively). The expected range of val¬ 
ues found for the models of Type-1 AGN is almost identical to 
the range of values found for PGQSOs and Sis. Moreover, the 
steepness of the spectra in the model of Type-2 AGN is expected 
to include the same range of values than Sis but extended toward 
lower values. Thus, Type-2 AGN are expected to show steeper 
spectra than Type-1 AGN. This is in agreement with our re¬ 
sults, where S2s tend to show lower log(yLy(20yum)/yLy(30yum)) 
than Sis. Seven out of the 20 SBs have a steepness fully 
consistent with the model of Type-2 AGN. All the objects 
below the minimum steepness predicted for Type-2 AGN 
by the Clumpy libraries can be considered as SB-dominated 
(log(yLy(20yum)/yLy(30yum)) = -0.24). However, there are SBs 
above that limit. Thus, log(yLy(20yum)/yLy(30yum)) < -0.24 in¬ 
dicates that the spectrum is SB-dominated but we cannot discard 
that spectra showing log(yLy(20yum)/yLy(30yum)) > -0.24 might 
also be SB-dominated. 

We have classified as SB-dominated spectra those showing 
log(yLy(20yum)/yLy(30yum)) < -0.24. Among the S2s, six ob¬ 
jects are therefore SB-dominated, three of them already classi¬ 
fied as SB-dominated according to the strengths of the PAHs 
and silicate features. Combining both methods together, ten out 
of the 31 S2s (32%) are SB-dominated. 

Thirteen out of the 41 LINERs show 
log(yLy(20/im)/yLy(30yum)) < -0.24. Two are LINERls 
and the remaining 11 are LINER2s. Among them, eight were 
already classified as SB-dominated using the strengths of PAHs 
and silicate features. Interestingly, only two LINERs classified 
as SB-dominated with the diagram seen in Eig. are not 
SB-dominated using the steepness of the spectra. Eifteen out of 
the 40 LINERs (37.5%) show signatures of being SB-dominated 
once the two methods presented in this section are considered 
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Fig. 5. (Left): average spectra for PGQSOs (orange), Sis (red), AGN-dominated S2s (yellow), and AGN-dominated LINERs (purple). (Right): 
average spectra for SBs (green), SB-dominated S2s (yellow), and SB-dominated LINERs (purple). We also show one standard deviation as a 
shaded region using the same colours. Each spectrum is also marked with different symbols at 21pm for clarity of the plot: PG QSOs (orange star). 
Sis (red square), SBs (green triangle), S2s (yellow square), and LINERs (purple circle). The average spectra are scaled to the flux at 15pm. 


together. The fraction of SB-dominated LINERs is similar to 
that of S2s. 


4.3. Goodness of the methodology to trace nuclear 
properties 

We have considered a mid-IR spectrum as AGN-dominated 
if it obeys these three criteria: EW(PAH 6.2pm) < 0.233yum, 
T 9.7 < 1.25, and log(yLy(20yum)/yLy(30yum)) > -0.24. To study 
the goodness of this method to select AGN-dominated 
Spitzer/IRS spectra, we have plotted in Eig. the \2pm flux 
for the SpitzerjlRS spectra versus the same quantity for ground- 
based measurements for the 63 objects for which these mea¬ 
surements are available (see Section |^. Ground-based and 
SpitzerjfRS \2pm fluxes show a linear relation although the 
dispersion is high (Pearson correlation coefficient of r=0.64, 
P(null) = 1.6 X 10“^). Moreover, the slope of the best linear fit 
(dashed fine in the left panel of Eig.|^ is flatter than the one-to- 
one relation (continuous line). 

Most of the SpitzerjfRS spectra with larger 12yum flux than 
those from ground-based measurements are SB-dominated ac¬ 
cording to the EW(PAH 6.2pm) method (see middle panel of 
Fig. 0. However, not all the SB-dominated SpitzerflRS ob¬ 
jects show a \2pm flux excess in the SpitzerflRS spectra com¬ 
pared to the ground-based measurements. We rule out the ex¬ 
planation of a distance effect in which more distant objects 
might include more SB-contribution in the nuclear spectra, be¬ 
cause all our objects are nearby and no particular trend is found 
comparing SB- and AGN-dominated sources. Alternatively, this 
result might have two explanations: ( 1 ) our method to select 
SB-dominated spectra is too restrictive, and could include SB- 
dominated spectra that are actually AGN-dominated at \2pm, 
and ( 2 ) the SpitzerflRS spectra do not contain extra emission 
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compared to the ground-based measurements and both are trac¬ 
ing a nuclear SB-dominated spectrum. 

If we select only AGN-dominated Spitzer spectra, the cor¬ 
relation between these two quantities improves (see Eig. 
right panel) with a Pearson correlation coefficient of r=0.88 
(P(null) = 5.7 X 10“^"^). The linear fit to the AGN-dominated 
sources (dashed fine in Eig. right panel) is very close to the 
one-to-one relation (continuous fine in Eig. right panel). The 
only two outliers are NGC 4594 and NGC 5866. Thus, when the 
AGN-dominated spectra are selected the nuclear 12pm flux ob¬ 
tained with ground-based data is very close to the value obtained 
by the Spitzer/IRS spectra. This reinforces our methodology as 
a good tool to isolate AGN-dominated mid-IR Spitzer/IRS spec¬ 
tra. 

5. Average spectra 

Eig. shows the average spectrum for each class of objects. 
These average spectra have been computed after normalising 
them to the flux at 15pm. The shaded regions show the stan¬ 
dard deviation over the average spectrum. We have computed 
the mean value for S2s and LINERs according to our mid-IR 
classification (see previous section) as AGN-dominated and SB- 
dominated (left and right panels, respectively). 

The average spectra of Sis and PGQSOs show similar 
shapes, showing the silicate feature in emission and similar 
steepness of the spectra (see left panel in Eig. [^. The relative 
differences between these two classes are an enhancement of the 
strength of the silicate feature in emission for PGQSOs com¬ 
pared to Sis and a slightly steeper spectrum for Sis compared 
to PG QSOs. The average SB spectrum is very different to that 
of Sis or PGQSOs (see right panel of Eig. [^. The main dif¬ 
ferences are strong PAH features, a steep spectrum, and deep 
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Fig. 6. The average spectra of AGN-dominated LINERs according to 
different subclassifications. Erom top to bottom: (a) LINERls (light ma¬ 
genta) and LINER2s (dark magenta); (b) bright LINERs (dark grey) and 
faint LINERs (light grey); (c) Compton-thin (light green) and Compton- 
thick candidates (dark green); (d) objects classified at X-rays as AGN 
(dark blue) and non-AGN candidates (cyan). The average spectra for 
Sis (red) and AGN-dominated S2s (yellow) are also shown for compar¬ 
ison purposes. 


silicate absorption features. Moreover, the classical lines associ¬ 
ated with AGN emission such as [Ne V] at 14.3yum and 24.3/im 
or [O IV] at 25.9/im, are clearly detected in the average spectra 
of PGQSOs and Sis but are absent in the average SB spectrum. 
Note here that the association of these lines with AGN emission 
have been questioned, findin g them in some SB galaxies (see 
[Pereira-Santaella et al.|2Q10| ). 

The average SB-dominated spectrum for S2s and LINERs 
(right panel of Fig.[^ are very similar to that of the SBs. The av¬ 
erage AGN-dominated (left panel in Fig. and SB-dominated 
(right panel in Fig. spectra for S2s and LINERs are clearly 


different. This support our method as a good diagnosti c of SB- 
dominated mid-IR spectra. This was also suggested by |Alonso- 


Herrero et al.| ( |2014| ), finding that the Spitzer spectra of Sis and 


S2s are similar only if spectra with deep absorption silicate fea¬ 
tures are excluded from the analysis. The [Ne V] at 143pm 
and [O IV] at 25.9yum emission lines are clearly seen in the av¬ 
erage spectra o f SB-dominated S2 s and LINERs. This was al¬ 
ready shown by pudik et al.| ( [2009| ), finding these lines in a large 
fraction of LINERs. Thus, the average SB-dominated spectra of 
LINERs and S2s might still show mid-IR signatures of AGN na¬ 
ture, although the overall mid-IR spectra is not dominated by the 
AGN. 


The average S2 AGN-dominated spectrum (yellow spectrum 
in the left panel of Fig.|^ does not mimic Sis or PG QSOs. This 
average spectrum is steeper than those of S Is or PG QSOs. It also 
shows the silicate features in absorption while Sis and PG QSOs 
show an average spectra with silicate features in emission. This 
is expected since the silicate feature at 9.7pm and ISyum are 
predicted to be i n emission for Type-1 AGN and in absorption 
for Type-2 AGN ( [Nenkova et al.|200^ . These 
already been confirmed by observations (e.g. 

On the similarities, the average (AGN-dominated) S2 spectrum 
shows the presence of [Ne V] at 14.3 and 243pm and [O IV] at 
25.9yum emission lines. 

The average AGN-dominated LINER spectrum (purple spec¬ 
trum in the left panel of Fig. can be clearly distinguished from 
PGQSOs, Sis, S2s, and SBs. The main characteristic of this 
average spectrum is the rather flat continuum all over the 6-28 
pm wavelength range. Moreover, it shows strong PAH features 
at 11.3yum and 11pm. The [O IV] at 25.9yum emission line is 
prominent as in Sis, PGQSOs, and S2s. However, the [Ne V] 
at 14.3yum emission line is clearly undetected as is the [Ne V] 
at 24.3yum. Below 20pm this spectrum resembles that of SBs. 
However, it can be clearly distinguished from SBs because the 
average spectrum of LINERs do not show a steep spectrum and 
it lacks of the silicate absorption features seen in SBs. More¬ 
over, LINERs also show the [O IV] at 25.9yum that the SBs do 
not show. 


predictions have 
IShi et al.||2006|). 


5.1. Sub-classes of LINERs 


As explained in the introduction, the LINERs are a heterogenous 
family of objects. In order to study the subclasses of LINERs, 
Fig.j^shows the average spectra for several subclassifications of 
LINERs. Note that these average spectra have been computed 
including only AGN-dominated spectra. 

The average spectrum of objects optically classified as 
LINERls (panel (a), light magenta spectrum in Fig.[^ is steeper 
than that of LINER2s (panel (a), dark magenta spectrum). In 
fact, the average spectrum of Type-1 LINERs is consistent with 
that of AGN-dominated S2s (yellow spectrum). However, the 
dispersion in the average spectrum of LINERls is quite large 
because, among the five AGN-dominated LINERs 1, NGC4450 
shows a flat spectrum which is not consistent with the average 
trend for this class. 

We have also classified the AGN dominated LINERs into 
two classes attending to their Lx(2 - 10 keV): bright LINERs for 
objects with Lx(2 - 10 keV) > lO'^^erg s“^ and faint LINERs for 
those with X-ray luminosities below that limit. This is the limit 
where the torus is expect ed to disappear at a bolometr ic luminos¬ 


ity of Lboi ~ lO'^^erg s ^( Elitzur & Shlosman 2006| p^ Among 


This assumes a conversion between the X-ray lu minosity an d the 
bolometric luminosity of Lboi - 10xLx(2- 10 keV) ( |Ho|2008| ). This 
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Fig. 7. The 2-10 keV luminosity versus the I2pm luminosity, both in logarithmic scale for PGQSOs, Sis and S2s (a), SBs (b), AGN-dominated 
LINERs (c) and S B-dominated LINER s (d). The continuous-red and dashed-green lines show the best fit correlations for AGN and SBs, respec¬ 
tively, reported by Asmus et"n^p01l| ). The typical error for these measurements is shown as a cross in the top-left corner of panel (c). Errors for 
the X-ray luminosity are estimated as 10% of its value. The dot-dashed blue line and long-dashed light green line show the linear fit for PG QSOs, 
Sis, and S2s and for SBs, respectively. Grey arrows mark objects with reported upper-limits on the X-ray luminosity. 


the AGN-dominated LINERs, 7 are classified as bright LINERs 
and 18 are classified as faint LINERs. The resulting average 
spectra of bright and faint LINERs (dark and light grey spec¬ 
tra in panel (b) of Eig. are quite similar to that of LINER Is 
and LINER2s. Thus, bright LINERs show a steeper spectrum 
compared to faint LINERs, compatible with Sis. Note that the 


conversion factor is the most conservative value we have found (|Ho| 
[2009] ). Note that any higher value could include more objects as bright 
LINERs (see also the discussion). 
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dispersion of the average spectrum of bright LINERs is much 
lower than that of LINERls. We have investigated if a differ¬ 
ent morphology of the host galaxy for bright and faint LINERs 
could produce the differences in their average mid-infrared spec¬ 
tra. Interestingly, NGC 4450 is the only AGN-dominated Type-1 
LINER hosted in a late-type galaxy. Only four (out of the 18) 
faint LINERs and 1 (out of the 7) bright LINERs are hosted in 
late type galaxies (i.e. t>l). This lack of late type galaxies host¬ 
ing LINERs is expected since they are usually hosted in early 
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Fig. 8. The 2-10 keV luminosity versus the 12/im luminosity, both in logarithmic scales for Compton-thin (left), Compton-thick (middle), and 
the full sample after Compton-thickness correction (right). The typical error for these measurements is shown as a cross in the top-lef t comer of 
the right p anel. The continuous-red and dashed-green lines show the best fit correlations for AGN and SBs, respectively, reported by |Asmus et| 
[aL] ( |201l| ). The long-dashed blue line shows the linear fit to the full sample. Grey arrows mark objects with reported upper-limits on the X-ray 
luminosity. Small double arrows in the middle panel show the expected locus if the real intrinsic 2-10 keV luminosity is 10-70 times higher than 
that estimated. 


type galaxies ( [Carrillo et al.|1999| ). Thus, we have not found any 
particular tendency in their host-galaxy morphology. Indeed, the 
average spectrum of faint LINERs hosted in late type galaxies 
is consistent with the average spectrum of faint LINERs (and 
with the average spectrum of faint LINERs hosted in early type 
galaxies). 

Dark and light green spectra shown in panel (c) of 
Eig. show the average spectra for Compton-thin and 
Compto n-thick candidates reported by [Gonzalez-Martm et al.j 
( |2009B| ). Compton-thick sources are expected to be so ob¬ 
scured (Nh > 1.5 X lO^^cm"^) that the bulk of the X-ray con¬ 
tinuum emission of the AGN is only seen above 10 keV. Sources 
were classified according to the ratio between the fluxes of 
[O III]T5007A and the 2-10 keV fluxes, the EW of the EeKa 
emission line, and the hard X-ray slope of the spectrum. Six¬ 
teen of them were classifled as Comp ton-thin LINERs and nine 
were classifled as Compton-thick candidates. We do not And 
signiflcant differences between the average mid-IR spectra of 
Compton-thin and Compton-thick candidates. 

[Gonzalez-Martm et al.[ ( [2009 A[ ) classifled 22 out of the 25 
AGN-dominated LINERs at X-rays: 16 AGN candidates and 6 
non-AGN candidates. Both average spectra are quite flat (cyan 
and blue spectra in panel (d) of Eig.[^. Interestingly, PAH fea¬ 
tures and [Ne II] at 12.8yum strengths are larger for non-AGN 
candidates than for AGN candidates suggesting the dominance 
of the host galaxy contribution. 


6. X-ray versus mid-IR luminosities 

In AGN, the X-ray emission is thought to originate in the inner¬ 
most regions of the accretion flow by Comptonisation processes. 
The mid-IR emission is generally dominated by thermal emis¬ 
sion by dust at parsec-scale distances from the SMBH. The tight 
correlation found between the X-ray and mid-IR emission of 
AGN supports this physical connection, regardless of the Seyfert 
type ([Krabbe et al.][200It [Gandhi et al.[[2009| [Levenson et al.[ 
[2009[ ). Eurthermore, SB nuclei also show a relation between X- 
ray and mid-IR luminosities but it is offset when compared to 


that of AGN ( [Asmus et al.|20Il| . The study of this correlation 
could bring light on the dominant process involved in LINERs. 

Eig. [T] shows Lx (2 - 10 keV) versus yLy(I2yum) for our sam¬ 
ple. Note that yLy(I2yum) is computed using the Spitzer spectra. 
Sis, S2s, and PG QSOs (panel (a)) follow the same relation pre¬ 
viously found for AGN (continuous red line, from [Asmus et al.[ 
[2011[ ). The Pearson correlation coefficient for them is r = 0.82 
(P(null) = 4.7 X 10“^^) and the best linear fit (blue dot-d ashed 
line) is very close to that reported by [Asmus et al.[ ( [MTT] ) using 
high spatial resolution mid-IR images. Therefore, this correla¬ 
tion holds for most Sis, S2s, and PGQSOs even with the rela¬ 
tively low spatial resolution of the Spitzer spectra. This may 
suggest that for these sources the AGN continuum dominates the 
Spitzer[[RS emission. Only two S2s (NGC 3621 and NGC 4725) 
fall close to the correlation found for SB nuclei (dashed green 
line, also from [Asmus et al.|2Qll[ ). In our study NGC 3621 was 
classified as SB-dominated and NGC 4725 as AGN-dominated 
(see Section]^. SB nuclei (green up-side down triangles in the 
panel (b) in Pig.|7]) also fall into the expected correlation for them 
with a correlation coefficient r = 0.83 (P(null) = 4.5 x 10“^^). 
None of them are compatible with the AGN correlation. Thus, 
this correlation seems to be very effective for distinguishing pure 
AGN emission from SB emission. 


The bottom panels of Pig. [7] show Lx(2 - lOkeV) ver¬ 
sus yLy(I2yum) for LINERs. We split the plot into the AGN- 
dominated (panel (c)) and the SB-dominated (panel (d)) LIN¬ 
ERs. Most LINER nuclei are placed between the AGN and the 
SB linear relations. The linear correlation for them is not signifi¬ 
cant (r = 0.68 and P(null) = 8 x 10“^). Low luminosity LINERs 
are near SB correlation and high luminosity LINERs are parallel 
to the AGN correlation. This could be due to an underestimation 
of the X-ray luminosity (see below). 

Some LINERs behave as AGN while some others are similar 
to SBs (see Section]^. Many of the LINER nuclei classified as 
SB-dominated are placed along the SB correlation (all the Type- 
1 LINERs). Again, this confirms our method described in Sec¬ 
tion]^ as a good tool to discriminate SB- from AGN-dominated 
spectra. However, three LINER2s (NGC 3079, NGC 6240, and 
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NGC7130) are placed in the AGN correlation but were classi¬ 
fied as SB-dominated in the mid-IR. Interestingly, all of them are 
known Compton-thick AGN. This suggests that the 12/im fiux 
could be AGN-dominated although the full Spitzer/IRS spectrum 
is SB-dominated. Therefore, irrespective of the mid-IR spectral 
shape, the I2pm fiux is a good tracer of the AGN power (e.g. 
|Gonzalez-Martm et al.|2013| ). 

The inclusion of LINERs in this correlation allows to vali¬ 
date it at lower luminosities (Lx(2 - 10 keV) < lO'^^erg s“\ see 
also [Mason et al.||2012| ). The final linear correlation for AGN- 
dominated spectra is: 


log(Lx) = (-12.34 + 0.05) -r (1.26 + 0.01)log(yLy(l2pm)) (2) 

which is very significant (r = 0.92 and P(null) = 2x 10“^^). 
However, the slope of this correl ation is steeper (~ 1.26) than 
that previously found (~1.06, see Asmus et al.||2011| ). This ex¬ 
cess in mid-IR luminosity for faint LINERs was already found 
by [Mason et al.| ( [20T^ . They argued that this discrepancy could 
be due to optically thin material that obscures the inner parts of 
the AGN because most of them showed silicate features in emis¬ 
sion. However, our faint LINERs do not show in average silicate 
features in emission (see Section [5T] and Eig.[^. 

The most natural explanation is that the Compton-thick 
nature of some low-luminosity AGN results in an underes¬ 
timation of the true Lx(2 - 10 keV) of these sources. These 
Lx (2 - 10 keV) estimates come from studies using the spec¬ 
tra of LINERs at energies below 10 keV. However, Compton- 
thick sources show the bulk of the AGN power at energies 
above 10 keV. The intrinsic luminosity could be 10-70 times 
higher than the estimated using only energies below 10 keV in 
the Compton-thick scenario ( [Maiolino et al.[p^98) ). [Gonzalez^ 
[Martin et al.[ ( [2009B[ ) classified around 50% of their LINER 
sample as Compton-thick candidates. Lig. shows the AGN- 
dominated objects in our study attending to their Compton-thin 
(left) and Compton-thick (middle) classification. Most of the 
Compton-thin sources are close to the previously reported cor¬ 
relation for AGN. Compton-thick S2s are nicely placed along 
the AGN correlation found by [Asmus et al.[ ( [2011[ ). This is 
expected because the X-ray luminosities included for S2s are 
all corrected for their Compton-thick nature (most of them in¬ 
cluded in [Goulding et al.[[MT^ , either using X-ray measure¬ 
ments above 10 keV or assuming a factor between the observed 
and intrinsic X-ray luminosity for other Compton-thick AGN 
(see [Panessa et al.[[2006| for details in this conversion factor). 
Compton-thick LINER candidates tend to be shifted toward X- 
ray luminosities lower than predicted for the AGN correlation. 
Note, however, that most of them are not consistent with the SB 
correlation either. The double arrows of Lig. (middle panel) 
show their expected locus if the X-ray luminosity were ~10- 
70 times higher. Most of the Compton-thick LINERs can be 
placed in the AGN correlation if the correction is applied. The 
linear fit to the entire sample (long-dashed blue line in Lig. 
[^, once the intrinsic X-ray luminosity is corrected results in 
Lx(2 - 10 keV)(intrinsic) = 40 x Lx(2 - 10 keV)(observed) for 
Compton-thick LINERs (marked as white stars in Lig. right 
panel), which is very close to the linear relation found for AGN. 
Nonetheless, four LINERs and one S2 still remain very close to 
the SB correlation. 


7. Discussion 

The nature of LINER nuclei has been extensively studied since 
they were firstly discovered by Heckman] ( [1980[ ). Using multi¬ 
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wavelength information and several techniques as compactness, 
hardness, and variability, we now know that a large fraction of 
them host an AGN ([Maoz et al.[[2005[ [Gonzalez-Martm et al.[ 
[2009A|Bt[H5p008]). However, what makes them a unique class 
is still unknown. Singh et ah] ( |2013[ ) showed that post-main se¬ 
quence stars might be an important contributor to the optical fre¬ 
quencies. Interestingly, they argue that this population might be 
present in more powerful AGN, although it is outshined by the 
AGN itself. At X-rays, several authors have pointed out to ob¬ 
scuration as one of the main ingredients for their different na- 
ture( [Dudik et al.[|20Q9l [Gonzalez-Martm et al.[[2009B| ). There¬ 
fore, mid-IR frequencies are key to study this obscuration since 
the emission absorbed at optical and UV frequencies is expected 
to be reprocessed at those wavelengths. 

We have compiled a sample of 40 mid-IR spectra of LINERs 
observed with Spitzer/IRS and have compared them to samples 
of SBs, S2s, Sis, and PGQSOs. Although the low spatial res¬ 
olution of these Spitzer/IRS spectra is a disadvantage, we have 
been able to isolate SB-dominated from those which are not by 
using well known mid-IR tools (see Sect.]^. We focus this dis¬ 
cussion in two main issues about the obscuration of LINERs in 
light of the present results: (1) Compton-thickness and (2) torus 
signatures in LINERs. 


7.1. Compton-thickness 


Gonzalez-Martm et al.[ ( [2009^ showed that a large fraction (up 


to 50%) of their LINER sample might be Compton-thick, i.e. 
with such a high obscuration that the intrinsic continuum of the 
AGN is fully suppressed at energies below 10 keV. This was 
done using indirect arguments as the [O III] 45007A and to X- 
ray fiux-ratio or the EW of the neutral LeKtr line at 6.4 keV. 
However, this is not yet confirmed with direct observations be¬ 
cause LINERs are too faint to be observed above 10 keV with 
the open X-ray instrumentation to the communit}p] 

We have take advantage of using the X-ray to mid-IR cor¬ 
relation found for AGN to study the Compton-thick nature of 
LINERs. LINERs classified as Compton-thick candidates by 
[Gonzalez-Martm et al.[ p009B[ ) are systematically located be¬ 
low the relation found for AGN (see Lig.[^ middle panel). This 
is naturally explained if they are indeed Compton-thick AGN, 
since their X-ray intrinsic luminosity is underestimated by a fac¬ 
tor of 10 or more. This was already found by [Mason et al.[ ( |20T^ 
for a small number of objects. Lour additional LINERs and one 
Seyfert also have lower X-ray luminosities that expected accord¬ 
ing to their mid-IR luminosities (see Lig. left panel). Thus, 
either they host a Compton-thick nucleus or are SB-dominated. 
This suggests that the actual fraction of Compton-thick LINERs 
may be even higher than previously inferred ( [Gonzalez-Martm[ 
[et al.[[2009B[ ). This large fraction of obscuration might explain 
why post-main sequence stars dominate the optical spectrum in 
LINERs ( [Singh et~ar|[2013[ ), while the AGN completely domi¬ 
nates the emission in powerful and/or less obscured AGN. 

Observationally, a large fraction of A GN in the local Uni- 
verse are obscured by Compton-thick gas ( [Maiolino et al.|19^ 
[Matt et al.|2000[ ). Lrom the theoretical point of view, a sizeable 
population of mildly Compton-thick sources is postulated in all 
the AGN synthesis models for the X-ray background in order 
to match the intensity peak of the XRB spectrum at about 30 
keV ( [Comastri|[2004[ ). These theoretical and observational evi¬ 
dence are consistent with the results of the present paper. Lur- 


NuSTAR satellite is able to observe these faint source but it is part of 
a close collaboration. 









































































O. Gonzalez-Martm et al.: On the LINER nuclear obscuration, Compton-thickness and the existence of the dusty torus 



Fig. 9. Best-fit to clumpy models (dot-dashed light-blue line) for 
the average spectra, from top to bottom, of Sis (panel (a)), S2s 
(panel (b)), bright LINERs (panel (c)), and faint LINERs (panel (d)). 
Note that bright LINERs are those with an X-ray luminosity above 
Lx > 1 X lO'^^erg/s. All the spectra are normalised to their flux at I5jum. 
The shaded light-blue area shows the lower and upper bounds obtained 
with the clumpy models (see text). 


thermore, the fraction of Compton-thick Seyferts is lower than 
that reported in LINERs (and confirmed in this analysis). This is 
fully consistent with the paradigm in which the fraction of ob¬ 
scured sources increases when the luminosities decrease in the 
high-r edshift Universe (|Ueda et al.||2001| ) and in the local Uni¬ 
verse ( [Shinozaki et al.|2006| ). 


mean in LINER Is are large indicating that probably they are 
not a well defined class. The best distinction between two sub¬ 
classes is found for bright and faint LINERs, above and below 
Lx > 1 X lO^^erg/s, respectively (Eig.[^. 

The mid-IR spectrum of AGN is expected to be dominated 
by dust re-emission of optical/UV emission. Under the unified 
model of AGN this dust is distributed in a dusty torus ( |An-| 
|tonucci||l993| ). This torus was first postulated to be a smooth 
doughnut-like structure. However, nowadays the preferred sce¬ 
nario is a clumpy structure, i.e. made of dusty clouds within a 
toroid a l distribution ([Krolik & Begelm^|1988[ [Pier & Krolik 


7.2. Torus in LINERs 

We have found that the average mid-IR spectrum of LINERs 
is markedly different from those of SBs, Sis, S2s, or PGQSOs 
(see Pig. 1^. The former shows a rather flat continuum from 6 
until 28yum. On top of this flat continuum, prominent emission 
lines of [Ne II], [Ne III], [S III], and [OIV] are clearly seen. The 
[O IV] emission line suggests the presence of AGN emission 
although it could be explained by another emission mechanisms 
(see [Pereira-Santaella et al.||201()| ). The PAH features are also 
clearly seen. They indicate a non-negligible contribution of star- 
formation to the mid-IR spectra of LINERs. However, note that 
star-formation cannot be the dominant process for them because 
the SB mid-IR continuum is much steepef^ (see right panel of 
Fig.|g. 

We have also found that the shape of the average mid-IR 
spectrum of LINERs does not depend on its Compton-thickness 
or AGN nature at X-rays (see Pig. It depends on the opti¬ 
cal type, i.e. the average mid-IR spectrum of LINERls is dif¬ 
ferent to that of LINER2s. However, the scatter around the 


Note that in this discussion SB-dominated LINERs are excluded, see 
Section 1^ 


|1992| ). [Nenkova et ai.| ( 2002| ) developed a formalism for han¬ 
dling radiative transfer in clumpy media and applied it to the IR 
emission from the AGN dusty torus. [Nenkova et al.| ( |2008| ) ap¬ 
plied this formalism to develop a set of models for this clumpy 
medium. They showed that clumpy torus models are consistent 
with current AGN observations if they contain a number of dusty 
clouds along the equatorial axis of No ~ 5 - 15, each with an op¬ 
tical depth of Tv ~ 30 - lOO. jRamos Almeida et al.| ( 2011|) mod¬ 
elle d the IR SED for a sample of Sis and S2s by Nenkova ^ 
[ST] ( |2008| ) dusty torus models, using a Bayesian inference tool 
called BayesClumpy ( [Asensio Ramos & Ramos Almeida|2009 ). 
They found that the IR SED of both Sis and S2s are well fitted 
by these clumpy torus models although their intrinsic properties 
are different (larger covering factor for S2s than for Sis, see also 
[Alonso-Herrero et al.|2011j ). 

We have used the BayesClumpy tool to investigate whether 
clumpy models could fit our average mid-IR spectra. The me¬ 
dian best fit models found for the average mid-IR spectra are 
shown in Pig. as dot-dashed light-blue line together with 68% 
confidence intervals for all the parameters as light-blue filled re¬ 
gions (see [Ramos Almeida et al.||2009 for more details on th e 
modelling). As already shown by [Ramos Almeida et al.[ ( [2QTT] ), 
the continuum shape of Sis and S2s are very well represented 
by clumpy models (panels (a) and (b)). Bright LINERs (i.e. 
Lx(2 - 10 keV) > 1 x lO'^^erg/s, panel (c)) can also be repro¬ 
duced by these models. However, faint LINERs (panel (d)) are 
poorly represented by them. This may suggest that the optically 
thick torus emission may not longer dominate the mid-IR emis¬ 
sion of faint LINERs. 

In the last years, several pieces of observational evidence 
have supported the scenario in which a single continuum dis¬ 
tribution of clouds within a wind is responsible for both the 
BLR and the dusty torus ( [Elitzur & Shlosm^[2006| and ref¬ 
erences therein). Under this scenario, the difference between 
the BLR and dusty torus is just a change on the composi¬ 
tion of this wind at the dust sublimation radius. An imme¬ 
diate consequence of the disk-wind scenario is the prediction 
that the torus and the BLR disappear at bolometric lum inosi- 
ties below Lboi ~ lO'^^erg/s ( [Elitzur & Shlosm"^ 20Q6| . This 
limit is consistent with our results assuming a conversion fac¬ 
tor between the bolometric luminosity and X-ray luminosity of 
Lboi/Lx(2 - 10 keV) ~ 10 (see [Ho|2008| ). 

7.3. Alternative source of emission at mid-IR 

If faint LINERs are no longer dominated by the torus emission, 
what is the mechanism responsible for the mid-IR emission? 
This is a question for which we do not have a clear answer. Plau¬ 
sible contributors to the mid-IR emission are: (1) jet emission in 
the form of synchrotron radiation; (2) host galaxy contributors; 
and (3) advection-dominated accretion flows (ADAPs). Note that 
we have excluded the dust shell models because we do not see 
emission silicate features in the average mid-IR spectrum of faint 
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Fig. 10. (Top): Mid-IR spectrum of faint LINERs (grey continuous line), SBs (green dotted line), elliptical galaxy NGC 1407 (cyan dot-dashed 
line), synchrotron emission of a jet (pink dot-dashed line), ADAE (red dashed line), and . (Bottom): Mid-IR spectrum of faint LINERs (grey 
continuous line), SBs (green dotted line), AGB star SMP SMC 018 (long-dashed blue line), post-AGB star IRAS 01005-^7910 (dot-dashed orange 
line), and PN SMPLMC 79 (short-dashed red line). 


LINERs, although it might be relevant for some of them (e.g. 
NGC 3998, [Mason et al. |2013] ). Fig. \T0\ compares the shape of 
the mid-IR spectrum of faint LINERs with several mechanisms 
discussed along this section. Note that all of them are scaled to 
the emission at 15yum. 


7.3.1. Host galaxy contribution 

We discard star-formation as the dominant emission since the 
SB-dominated LINERs were not considered in the analysis, al¬ 
though it can still be present in a fraction as suggested by the 
PAH emission. 

Alternatively, it could be linked to the diffuse emission com¬ 
ing from the galaxy. LINERs ar e generally found in elliptical 
galaxies (see |Carrillo et al.|1999| ). Indeed in very few cases faint 
LINERs are hosted in late type galaxies (see Section |5.1| ). To 
compare the mid-IR spectra of LINERs with elliptical galaxies 
we have chosen the SpitzerJRS spectrum of the prototypical el¬ 
liptical galaxy NGC 140'^ Although the general slope of faint 


We have chosen NGC 1407 because it is one of the few elliptical 
galaxies with full coverage of the mid-IR spectrum. 
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LINERs resembles that of an elliptical galaxy, it lacks of PAH 
features and emission lines seen in faint LINERs. 


[Singh et aT] ( |2013| ) found that a large fraction of the optical 
spectrum of many LINERs could come from post- asymptotic 
giant branch (post-AGB) stars. In order to deeply investigate this, 
Fig. [T^ (bottom) compares the spectrum of faint LINERs and 
SBs with different stages of the AGB stars. 

Our average LINER faint spectrum shows PAH features (7.7, 
8.6, 11.3 and 17 pm) together with emission lines (e.g. [Ne II], 
[Ne III], [S III], and [O IV]). This is very different for AGB 
stars (see the spectrum of an AGB star SMP SMC 18 in Fig.p^ 
bottom) which show featureless continuum-dominated spectra 
( [Yang et al.|2004t[Sloan et al.|2006| ). Post-AGB stars (also called 
photo planetary nebulae, proto-PNe) show some nebular emis¬ 
sion lines and/or PAH features (see the spectrum of post-AGB 
star IRAS 01005-1-7910 in Fig.[^ bottom). However, they show 
a peak on the continuum emission above 25yum and they lack of 
the [O IV] line at 26pm seen in faint LINERs ([Cerrigone et al. 
[2009] ). 


The average spectrum of faint LINERs is more similar to 
evolved carbon-rich PNe, where the continuum is less promi¬ 
nent but PAH features and highly ionised nebular lines (e.g. [O 
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IV]) are still present ( [Stanghellini et al.|2007| ). Fig.(bottom) 
shows the comparison between the average spectrum for faint 
LINERs and the prototype carbon-rich PNe SMPLMC79. Al¬ 
though less prominent, the continuum of these carbon-rich PNe 
is still present. |Stanghellini et al.| ( |2Q12| ) showed that the dust 
continuum of carbon-rich PNe is well fitted with a black-body 
model with a temperatures ranging T -100-160 K. Moreover, at 
least 10^ PNe are needed if they are the only responsible for the 
mid-IR luminosity of faint LINERs, si nce PNe typically show IR 
luminosities of Lir -lx lO^'^erg s“^ ( Stanghellini et al.|2012| . 
Alternatively, a combination of elliptical galaxies (top panel of 
Fig. [T^ and carbon-rich PNe could also reproduce the mid-IR 
average spectrum of faint LINERs. Finally, we have checked the 
[OIJ/Hof ratio for our sub-sample of faint LINERs, and in eight 
out of the 22 of them the ratio is to o high (i.e. [QI]/Ha > 0.25 ) 
to be explained by post-AGB stars ( |Cid Fernandes et al.|2004| ). 
Thus, none of these stages of AGB stars can fully accomplish for 
the observed features of faint LINERs, although a combination 
of several of them can explain the mid-IR spectrum for some of 
these faint LINERs. 

Interestingly, a combination of carbon-rich post-AGB stars 
and PNe could explain the SB-dominated spectra. Moreover, the 
average spectrum of SBs SB-dominated S2s and SB-dominated 
LINERs show very similar spectra (see Fig.[^. Thus, in all these 
cases a combination of carbon-rich post-AGB and evolved PNe 
might be relevant at mid-IR. This is in agreement with the find¬ 
ings by [Singh et aI.| ( [20T^ , where post-AGB stars dominate the 
optical spectrum when the AGN is faint. 


7.3.2. Advection-dominated accretion flows 


ADAF mechanisms have been largely claimed as the main re¬ 
sponsible for the SED of LINERs and LLAGN in general (e.g. 
[Nemmen et al.||MT4| and references therein). In these models 
the accretion disc is truncated at the inner parts and replaced by 
a hot, geom etrically thick, optically thin accretion flow ( jNarayanj 
et al.|1998] ). This flow is radiatively inefficient accreting material 


to the inner parts. ADAF models are quite complex in showing 
many observational parameters. In order to compare our average 
mid-IR spectra with these models we have chosen the ADAF 
best-fit to NGC 1097 reported by jNemmen et al.|(2Q14| ). ADAF 
model (see Fig. [T^ fails to reproduce the slope of the spectrum 
and also lacks of the PAH features and emission lines seen in 
faint LINERs. 


7.3.3. Jet emission 


LLAG N are generally radio loud according to their SED ( |Ho| 
[2008 ). Several authors have pointed out that the full SED 
could be explained as emission from the jet (e.g. NGC 1052, 
|Femandez-Ontiveros et al.||2012j ). [Mason et al.j ( [2013[ ) showed 
that jet emission is relevant for radio-loud LLAGN. To test this 
hypothesis, here we adopt a simplistic jet model based on in¬ 
ternal shocks adapted for AGN (see [Nemmen et al.|2006[[MT4 


and references therein). The parameters of the model are those 
used in [Nemmen et al.[ ( [2014[ ), assuming a power-law index dis¬ 
tribution with index p = 2. Note that we have assumed that the 
optically thin part of the jet is the one dominating the mid-IR 
emission. Jet emission (similarly to the ADAF models) shows an 
spectral index which is opposite to that observed in faint LIN¬ 
ERs (see top panel of Fig. [^. However, a combination of jet 
and other mechanisms could also explain the average spectrum 
of faint LINERs. 


In summary, the shape of the mid-IR continuum of faint LIN¬ 
ERs resembles that of elliptical galaxies, with a plausible contri¬ 
bution of carbon-rich PNe. However, it cannot be entirely de¬ 
scribed by any of the emission mechanism explained above. An 
AGN component might also be present at some level to account 
for lines like [O IV], which are sometimes suggested as indica¬ 
tive of AGN activity (see [Dudik et al.|2009] ). Note, however, that 
the emission mechanism producing the [OIV] line has been pro¬ 
posed to be unrelated to the AGN ( Pereira-Santaella et al.|2010 ). 
Indeed it is present in the spectrum of carbon-rich PNe (see the 
bottom panel of Fig. [TO. 


8. Conclusions 

We have analysed a sample of mid-IR spectra of 40 LINERs ob¬ 
served with SpitzerjXRS. We compared the LINER sample with 
PGQSOs, Seyferts (SI and S2), and Starburst (SB). The main 
results of this paper are: 

- We have developed a methodology to exclude SB-dominated 

mid-IR spectra based on the silicate optical depth the 

strength of the 6.2yum PAH feature, and the steepness of the 
mid-IR spectra. We have found that 25 out of the 40 LINERs 
do not show signatures of being SB-dominated. The fraction 
is similar to that obtained for Type-2 Seyferts. 

- LINERs fall into the X-ray to mid-IR correlation for AGN 
only if Compton-thick candidates proposed at X-rays show 
an X-ray luminosity around 70 times higher than that com¬ 
puted at X-rays. This is expected if they are confirmed 
as Compton-thick AGN (see Section [^. Furthermore, four 
LINERs not previously classified as Compton-thick candi¬ 
dates are also consistent with being Compton-thick accord¬ 
ing to the X-ray to mid-IR relation. 

- Bright LINERs (i.e. Lx(2 - 10 keV) > lO^^erg/s) 
tend to have an average mid-IR spectrum similar to 
that of Type-1 Seyferts. However, faint LINERs (i.e. 
Lx (2 - 10 keV) < lO^^erg/s) show a flatter average mid-IR 
spectrum, showing prominent emission lines. We suggest 
that this favours the disappearance of the dusty torus in LIN¬ 
ERs with bolometric luminosities below Lboi - lO'^^erg/s, 
as predicted theoretically. Therefore, LINERs with bolo¬ 
metric luminosities below Lboi - lO'^^erg/s might not longer 
be dominated by the torus in the mid-IR. Their mid-IR 
continuum emission resembles that of an elliptical galaxy 
although an AGN or a jet component together with some 
star-forming contribution (post-AGB stars and PNe) might 
also be present at some level. 


A final caveat is that a large fraction of LINERs seem to be 
Compton-thick at X-rays (see previous subsection) although the 
torus emission have disappeared for them. How is that possible? 
A natural explanation is that the Compton-thick material seen 
at X-rays is not produced in the dusty torus. Dusty material ab¬ 
sorbs continuum radiation both in the UV/optical and X-rays, 
and therefore the dusty torus might also be responsible for part 
of the X-ray obscuration. But dust-free gas attenuates just the 
X-ray continuum, so clouds inside the dust sublimation radius 
will provide additional obscuration only in this band. Conclusive 
evidence for such absorption comes from the short timescales 
for transit of X-ray absorbing clouds across the LOS, which es¬ 
tablish the existence of obscuri ng clouds inside the dust subli¬ 
mation radius (e.g. NGC 1365, [Risaliti et al.[[2009| ). A natural 
explanation is that the Compton-thick clouds seen at X-rays in 
faint LINERs are produced in dust-free regions within the dust 
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sublimation radius. Could LINERs be at a stage where the torus 
structure is already gone while the BLR is still present? This 
naturally explains why the optical spectrum of LINERs needs a 
population of post-main sequence s tars, blocking the source of 
ionization that produces NLR lines ( [Singh et al.|2013| ). Alterna¬ 
tively, this Compton-thick gas might not be related to the AGN, 
but rather to the host-galaxy as a result of galaxy interactions. 
In favour of that, many of the well known Compton-thick AGN 
are hosted in galaxy pairs or clusters (e.g. NGC 3690, |Goulding| 
jet al.|201^ [Gonzalez-Martm et al.|2013|) b eing even binary nu¬ 
clei (e.g. NGC 6240, [Komossa et al.|2ofe] ). A closer look at the 
environment of Compton-thick AGN needs to be undertaken in 
order to get more light on this possibility. 
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